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Abstract 


Steady  two-dimensional  stratified  shear  flows  of  a compressible 
fluid  over  a mountain  is  studied  with  the  long-wave  assumption 
(hydrostatic  balance) . For  constant  entropy  profile  far  upstream, 
we  obtained  a close  form  solution.  For  general  stratification  and 
shear,  the  speed  of  certain  points  in  the  flow  may  vanish.  We 
treated  these  points  as  the  birth-bed  of  flow  fields  with  closed 
streamlines  (recirculating  region) . Constants  of  motion  along  a 
streamline  are  used  to  construct  the  outer  boundary  of  these 
regions.  Numerical  solutiors  with  data  upstream  to  the  Rocky 
Mountains  were  constructed  with  the  recirculating  region  in  the 
flow  field. 


1.  Introduction 


The  flow  of  compressible  stratified  fluids  over  an  obstacle 
is  the  subject  of  this  paper.  Assuming  a steady  unperturbed 
state  far  upstream  of  the  obstacle  and  the  existence  of  either 
a free  or  rigid  surface  as  the  upper  boundary  condition  (see 
figure  1) , we  determine  the  resulting  flow  field  with  or  without 
closed  streanlines.  We  emphasize  the  introduction  of  closed 
streamlines.  It  is  shown  that  for  certain  upstream  conditions 
and  some  obstacle  height,  the  velocity  of  the  fluid  vanishes  at 
definite  points  within  the  flow  field.  It  is  at  these  points 
that  the  regions  of  closed  streamlines  originate.  For  simplicity 
we  restrict  ourselves  to  the  case  of  steady,  two  dimensional, 
isentropic,  compressible  flows  in  the  presence  of  gravity. 

In  section  2,  using  the  hydrostatic  or  long  wave  approxima- 
tion we  derive  a pair  of  first  order,  nonlinear,  ordinary 
differential  equations  which  govern  the  flow  in  this  problem. 

In  the  derivation  of  these  equations  we  obtain  three  constants 
of  motion  along  streamlines  which  prove  useful  in  the  analysis  of 
flows  with  closed  streamlines.  We  also  observe  that  the  assump- 
tion of  hydrostatic  balance,  along  with  the  condition  of  either 
a free  or  rigid  upper  boundary,  allows  us  to  solve  the  problem 
in  the  following  indirect  manner:  given  upstream  velocity  and 

entropy  profiles  and  the  velocity  of  the  fluid  along  the  upper 
boundary,  find  the  velocities  in  the  interior  of  the  flow  field 
and  the  height  of  the  obstacle.  In  fact,  the  flow  at  any  down- 
stream station  depends  only  on  the  conditions  far  upstream  of  the 
obstacle  and  on  the  height  of  the  obstacle  immediately  below. 


This  allows  us  to  construct  a master  solution  for  given  upstream 
conditions  which  contains  all  possible  steady  flow  fields  corre- 
sponding to  obstacles  of  height  less  than  or  equal  to  the  critical 
height,  which  is  the  maximum  obstacle  height  allowed  by  the  solu- 
tion. The  corresponding  procedure  for  incompressible  flow  without 
closed  streamlines  was  introduced  by  Lee  and  Su  (1977). 

In  section  3 we  assume  that  the  flow  is  homentropic  (i.e. 
that  entropy  is  constant  throughout  the  fluid)  and  show  that  in 
this  case  the  system  of  ordinary  differential  equations  derived 
in  section  2 can  be  solved  analytically.  For  a constant  velocity 
profile  far  upstream  of  the  obstacle,  we  obtain  a master  solution 
which  contains  all  possible  steady  flow  fields  corresponding  to  a 
particular  choice  of  the  Froude  number.  Yih  (1965)  constructed 
solutions  for  homentropic  flows  under  the  assumptions  that  specific 
energy  is  constant  throughout  the  fluid  and  that  the  dynamic 
compressibility  can  be  neglected.  Yih  (1965)  also  constructed 
solutions  for  flows  corresponding  to  slight  stratifications  in 
entropy  and  specific  energy  assuming  again  that  dynamic 
compressibility  may  be  neglected.  These  solutions  differ  from 
ours  in  that  Yih  did  not  make  the  long  wave  approximation. 
Consequently  his  method  yields  solutions  containing  lee  waves, 
while  our  method  docs  not. 

In  section  4 we  allow  variable  entropy  and  velocity  profiles 
far  upstream  of  the  obstacle.  In  this  case  analytic  solutions 
are  no  longer  possible  and  thus  we  resort  to  numerical  techniques 
for  the  integration  of  our  system  of  ordinary  differential 
equations.  Solutions  with  upstream  conditions  corresponding  to 
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flow  over  the  Uocky  Mountains  are  obtained.  These  solutions 
typically  contain  one  or  more  stagnation  points  in  the  interior 
of  the  flow  field.  Appropriate  use  of  the  three  constants  of 
motion  along  streamlines,  obtained  in  section  2,  provide  us 
with  three  constraining  conditions  to  construct  discontinuous 
solutions  containing  split  streamlines.  It  should  be  noted 
however,  that  our  analysis  does  not  provide  any  information 
concerning  the  flow  in  the  recirculating  region.  The  method 
used  to  study  the  splitting  streamlines  is  closely  related  to 
the  one  used  by  Margolis  and  Su  (197G)  in  their  analysis  of 
flows  with  a critical  layer. 


2 . Equations  of  Motion 

The  equations  governing  steady,  compressible,  isentropic, 
two  dimensional  flows  in  the  presence  of  gravity  are  the 
horizontal  and  vertical  momentum  equations 
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the  equation  of  continuity 

§x  (pu)  + 37  (pv)  = 0 ' 
and  the  equation  of  state 

p/pY  = A(s)  , 


(2.1) 

(2.2) 


(2.3) 


(2.4) 


where  s denotes  entropy  and  y is  the  ratio  of  specific  heats 


y=c  /c  . 
’ p v 


From  equations  (2.1) ,(2.2) ,(2.3)  and  (2.4)  it  can  be  shown 


that 
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where  <p  is  the  stream  function  defined  by 
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Using  the  hydrostatic  relation  we  replace  the  vertical 
momentum  equation  (2.2)  with 

3p 


Tfy 
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Equations  (2. 5),  (2. 6)  and  (2.7)  become 


A = A ( i|j) 

2 2 

H = H_+gy  + _£r  = H ( ) 

and  w = (j-  + gy)j|  + 337  = 
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(2.9) 

(2.10) 
(2.11) 


Following  Long  (1953)  we  define  a new  vertical  coordinate 
Y = Y ( tl» ) = Y (x,y)  as  the  height  far  upstream  of  the  streamline 
which  passes  through  the  point  (x,y)  (see  figure  2).  We  then 
have  that  the  horizontal  component  of  the  fluid  velocity  far 

^ I 

upstream  is  given  by  U(Y)  = — ~ , where  p (Y)  is  the  density 

pQ  3 Y o 

far  upstream.  Equation  (2.11)  now  takes  the  form 
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Introducing  t(Y)  = y-Y  we  have  that 

dC  ,Po,  U . 

dY  'p  ' u * 


(2.13) 


Thus  we  have  reduced  our  problem  to  that  of  solving  the  following 
pair  of  first  order  ordinary  differential  equations 
-2 
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where  the  non-dimensional  barred  variables  are  defined  as 

A = A/A  , U = U/U  , u = u/U  , y = y/H,  Y = Y/H,  p = p/p  and 

g g g 9 

pQ  = pQ/Pg.  Also  pQ  is  defined  as  the  pressure  far  upstream, 

H is  the  distance  between  the  upper  and  lower  surfaces  far  upstream, 

A , U and  p are  evaluated  far  upstream  at  Y=0  and  F,  the  Froude 
g g g 

2 

number,  is  given  by  F = U^/gH. 

Using  equation  (2.10)  we  derive  the  following  equation 
expressing  the  pressure  in  terms  of  velocity  and  the  displacement 


of  the  streamline  c(Y), 
y-1 


p(Y)  = [pQ  Y + 2~|  (U2-u2)  - Jc(Y)l 


_Y 
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(2.16) 


where  k = * In  genera-*-  WG  determine  pQ(Y)  using  the 

hydrostatic  relation  (2.8)  applied  to  the  flow  far  upstream. 


3. 


Analytic  Solutions  for  Homentropic  Flows. 


Assuming  that  the  flow  is  homentropic  (i.e.  constant  entropy 
throughout  the  flow  field),  the  solution  of  the  system  (2.14) 
and  (2.15)  can  be  represented  analytically.  The  hydrostatic 
relation  (2.8)  applied  to  the  flow  far  upstream  of  the  obstacle 
takes  the  form 


dY 


which,  due  to  the  assumption  that  A(Y)  = constant,  can  be 
integrated  to  yield 
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Equations  (2.14)  and  (2.15) 
(3.1)  and  the  assumption  of 
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become,  using  equations 
homentropy 


(3.1) 


(2.16)  and 


(3.2) 
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where  K = k + ~ F(U2-u2). 

The  solution  of  (3.2)-(3.3)  corresponding  to  a particular 
choice  of  U(Y)  and  satisfying  the  riqid  upper  boundary  condition 
y(l).=  1 is  given  by 
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where  we  have  used  u = 
of  (3.2).  The  value  of 
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by  the  obstacle  height  h immediately  below  this  position  through 
the  equation 

y (0)  = h 

and  a=0  corresponds  to  zero  obstacle  height.  The  allowable 
range  of  values  of  a is  restricted  b y the  requirement  that 
the  pressure  must  remain  positive.  Using  equation  (2.16)  we 
find  that  this  implies  that  a must  satisfy 

a < | (k-1)  (3.5) 

which  leads  to  the  further  restriction  that  k must  be  greater 
than  1 since  a=0  corresponds  to  unperturbed  flow. 

If  we  also  assume  that  U(Y)  is  constant,  the  solution  of 

(3.2)-(3.3)  for  the  case  of  a rigid  upper  boundary  becomes 

Y Y Y Y-l 

y(Y)  = E - {(E-l)y~1  + -J  [ (k-Y)  1-1  - (k-1 ) Y_1  ] } Y . (3.6) 

We  define  the  critical  Froude  number  F to  be  the  Froude  number 

c 

for  which  the  critical  height  is  zero.  It  is  easily  shown  that 
for  the  above  solution  (3.6) 
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(k-1) 
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(3.7) 


In  figures  3,  4 and  5 we  plot  the  obstacle  height  y(0)  versus  u 

for  the  three  qualitatively  distinct  cases  F < Fc , F = F^  and 

F > F . In  figure  6 we  plot  the  critical  height  h versus  F. 
c c 


.2. 


It  should  be  noted  that  there  is  no  added  difficulty  in  the 
case  that  the  upper  boundary  is  a free  surface.  The  appropriate 
boundary  condition  for  this  case  is  that  the  pressure  be  constant 
along  the  uppermost  streamline  Y=1  or  equivalently 

y(l)  = 1 + | F(U2-U2)  (3.8) 


where  we  have  used  equation  (2.16).  The  exact  solution  for  this 

case,  assuming  homentropy  and  constant  U(Y),  is 

Y Y Y Y~1 

y (Y)  = k - {(k-l)Y_1  + y.[(k-Y)Y-1  - (k-1 ) Y_1  ] } Y (3.9) 


and  the  critical  Froude  number  F is 
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(3.10) 


We  make  the  observation  that  our  system  of  ordinary  differential 
equations  (2.14)  and  (2.15)  approaches  the  equations  governing 
incompressible  flow  over  an  obstacle  as  y -+  °°  . Taking  this  limit 
we  note  that  Fc  -►  « for  the  case  of  a rigid  upper  boundary  and 
F -»■  1 for  the  case  of  a free  upper  boundary. 

The  above  results  are  easily  extended  to  the  case  where 
entropy  is  constant  in  two  or  more  horizontal  layers  which 
together  make  up  the  entire  flow  field. 


4.  Numerical  Solutions  for  Nonhomen tropic  Flows. 

For  arbitrary  velocity  and  entropy  profiles  far  upstream  of 
the  obstacle  we  are  no  longer  able  to  obtain  analytic  solutions 
Therefore  we  integrate  equations  (2.14)  and  (2.15)  numerically. 


The  procedure  that  was  followed  consisted  of  numerically 


integrating  the  initial  value  problem 
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4(1)  = 0 u(l)  = a 
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from  Y=1  to  Y=0  using  a 4th  order  Runge-Kutta-Gill  subroutine. 

The  parameter  a is  allowed  to  range  over  all  positive  numbers 
such  that  p > 0. 

Frequently  flows  of  this  type  contain  stagnation  points  in 
the  interior  of  the  flow  field.  In  order  to  perform  the  integra- 
tion of  equations  ( 2 . 14 ) - ( 2 . 15 ) at  positions  on  the  downstream 
side  of  the  stagnation  point  it  is  necessary  to  allow  the 
streamline  passing  through  the  stagnation  point  to  split  (see 
figure  7 ) . 

Thus  for  positions  on  the  downstream  side  of  the  stagnation 
point,  the  integration  is  done  in  two  stages.  First  we  integrate 
the  initial  value  problem  (2 . 14 ) , (2 . 15)  and  (4.1)  from  Y=1  to 
Y=Ys  where  Yg  denotes  the  height  far  upstream  of  the  streamline 
passing  through  the  stagnation  point.  In  order  to  complete  the 
second  stage  of  the  integration  (from  Y=Yg  to  Y=0)  we  derive 
three  jump  conditions  which  constrain  the  flow  in  the  region 
Y < Yg.  These  jump  conditions  follow  from  the  three  constants 
of  motion  along  streamlines  (2.9),  (2.10)  and  (2.11).  From  (2.9) 


v/e  have  that 


1 0 
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where  [(•)]  = (•)[  - ( • ) | • 
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Equation  (2.10),  along  with  the  hydrostatic  assumption  and  the 
assumption  that  entropy  is  constant  in  any  recirculating  region 
bounded  by  a split  streamline,  implies 

[u2]  = 0 . (4.3) 

The  assumption  of  constant  entropy  in  the  recirculating  region 
is  motivated  by  stability  considerations,  although  the  present 
analysis  cannot  provide  any  details  of  the  flow  in  this  region. 
Equation  (2.11)  and  the  above  results  (4.2)  and  (4.3)  imply 

f |y  Un  A1/y)  [y]  = [|y  (u2)]  (4.4) 

where  [y]  = y(Y+)  - y(Y  ) is  the  splitting  width  of  the  split 
streamline,  and  is  unspecified  in  this  analysis. 

This  procedure  was  carried  out  for  upstream  velocity  and 
entropy  profiles  corresponding  to  flow  over  the  Rocky  Mountains 
(see  figures  8 and  9).  Figure  10  is  a plot  of  velocity  u(l) 
along  the  uppermost  streamline  versus  the  obstacle  height  y(0). 
This  is  a representation  of  the  master  solution  for  this  problem. 
Figure  11  is  a plot  of  the  streamlines  for  the  flow  over  an 
obstacle  whose  height  is  equal  to  the  critical  height. 
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